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INTRODUCTION 

There  has  been  a  remarkable  effort  to  elucidate  the 
behavior  of  polymers  under  fatigue  loading  conditions. 

Recently,  the  phenomenology  of  polymer  fatigue  has  been 
thoroughly  reviewed  by  Hertzberg  and  Manson  [1].  In  analogy 
with  the  more  well-studied  low  molecular  weight  materials, 
fatigue  failure  in  polymers  occurs  due  to  the  accumulation 
of  damage  resulting  from  repeated  local  plastic  deformation. 

It  is  generally  understood  that  the  formation  of  crazes  and 
associated  cracks  is  responsible  for  the  ultimate  failure  of 
many  glassy  polymers.  However,  the  exact  nature  of  the 
mechanisms  involved  is  not  yet  well  understood. 

In  view  of  the  obvious  importance  that  affiliates  with 

the  prediction  of  Fatigue  Crack  Growth  (FCG)  rates,  various 

empirical  as  well  as  analytical  approaches  have  been  adopted 

in  an  attempt  to  develop  the  law  of  FCG.  Conventionally,  FCG 

"laws"  are  established  by  expressing  the  propagation  rate  of  a 

starter  crack  in  terms  of  the  applied  stress  or  a  related 

function.  It  is  not , therefore ,  surprising  to  find  several 

•  __  • 

such  empirical  "laws"  in  the  form:  £  =  f(K)  ,  where  l  is  the 
rate  of  crack  extension,  K  is  the  stress  intensity  factor,  and 
m  is  a  numerical  factor. 

As  demonstrated  by  numerous  examples  in  the  recent 
literature,  the  conventional  approach  suffers  from  the  lack 
of  applicability  to  experimental  data.  This  is  due  to  the 
fact  that  these  developments  rest  upon  the  classical  view  of 
fracture  mechanics.  That  is,  the  failure  process  involves  only 
the  propagation  of  a  single  crack-cut;  with  or  without  a  plastic 
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zone  ahead  of  its  tip.  However,  direct  observations  in  polymers 
[2,3]  and  other  materials  [4]  show  that  the  crack  propagates 
through  a  damage  zone.  Specifically,  Bevis  and  Hull  [3]  have 
shown  that  the  application  of  a  tensile  stress  to  a  precracked 
specimen  of  polystyrene  produced  a  halo  of  microcrazes 
surrounding  the  slowly  growing  root  crack  (Fig.  la).  A  similar 
zone  was  also  observed  under  low  cycle  fatigue  of  thin  samples 
of  polycarbonate  and  polystyrene  [5].  The  formation  of  a  zone 
of  microcrazes  around  the  crack  tip  could  be  a  source  of 
substantial  energy  dissipation.  This  can  be  considered  analogous 
to  the  classical  crack-tip  plastic  zone.  Clearly,  a  major 
alteration  in  the  law  of  crack  propagation  would  be  expected  due 
to  this  phenomenon. 

In  this  paper,  we  present  a  more  precise  analysis  of 
the  problem  of  fatigue  crack  propagation.  The  treatment  takes 
into  account  the  role  played  by  the  microcrazing  zone  and  its 
effect  on  the  fatigue  crack  behavior  under  sinusoidal  loading 
conditions.  A  more  general  expression  of  the  failure  process 
under  such  loading  conditions  is  formulated.  The  model  provides 
a  good  description  of  experimental  data  on  the  crack  growth  rate 
in  polystyrene  for  which  other  models  proved  inadequate. 

THEORETICAL  DEVELOPMENTS 

The  Crack-Crazing  Zone  (CCZ):  Based  on  the  observations 
outlined  in  the  introduction,  failure  processes  in  polymers, 
particularly  under  fatigue  loading  conditions,  involved  crack 
propagation,  root  craze  extension,  and  the  simultaneous  spread 
of  microcrazes  around  the  crack-root  craze  (Fig.  la).  The  "crack" 
system,  thus  described,  is  coined  as  "Crack-Craze  Zone"  and  is 


abbreviated  CCZ. 


Along  an  axis  perpendicular  to  the  CCZ  direction, 
microcrazing  density  assumes  a  bell-shaped  distribution. 
However,  for  simplicity,  the  microcrazing  distribution  is 


approximated  by  an  area  o-f  uniform  microcraze  distribution 
with  a  distinct  boundary.  The  assumed  boundary  separates  the 
microcrazed  zone  from  the  initial  material  (Fig.  lb).  We 
consider  CCZ  propagation  through  a  thin  polymer  sample  so  that 
two  dimensional  analysis  may  be  applied  and  that  isothermal 
conditions  may  be  considered.  Two  geometrical  parameters  of 
CCZ  are  accounted  for:  length  (£)  and  width  (a)  (Fig.  lb). 

Description  of  the  CCZ  history,  therefore,  requires  the 
definition  of  the  functions  £  =  £(t)  and  a  =  a(t),  where  t 
stands  for  time. 

Thermodynamics  of  CCZ  Propagation:  Viewed  in 
thermodynamic  terms,  the  propagation  of  CCZ  is  an  irreversible 
process  which  is  best  treated  by  methods  of  irreversible 
thermodynamics.  The  thermodynamics  potential  G(Gibbs  free 
energy)  for  the  entire  system  can  be  expressed  as  the  sum  of 
two  terms 

G  =  G1(o,  T,  £ ,  a)  +  G2(a,  T,  £,  a)  (1) 

G^  and  G2  are  the  thermodynamic  potentials  of  the  CCZ  and  the 
complementary  part  of  the  body.  The  stress  o,  and  temperature, 

T,  are  external  parameters,  whereas  £  and  a  are  internal 
parameters. 

Gibbs  free  energy,  on  the  other  hand,  may  be  resolved 
into  Helmholtz  free  energy  of  the  initial  material  Fq(T), 

Helmholtz  free  energy  of  deformation  F£  (o,T)  and  the  work  A 
done  by  external  forces.  Therefore,  we  may  write 

A 
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G(k)  Fo(k)(T)  +  Fe(k)(o,T)  ~  A(k)  (2) 

k  =  1,  2  indicates  the  CCZ  and  the  complementary  part  of  the 
body,  respectively.  Because  the  total  displacement  u  may  be 
considered  as  the  sum  of  a  reversible  part  (elastic)  u 

e 

and  irreversible  (plastic)  part  Up  which  are  associated  with 

corresponding  deformations  e  and  e  ,  the  work  A  may  be  expressed 

e  p 

as  the  sum  of  an  elastic  part  and  a  plastic  part,  i.e., 

A(k)  =  Ae(k)  +  Ap(k)  (3) 

By  definition,  the  potential  energy  (P)  associated  with  the 

elastic  deformation  is  given  by 

P(k )  =  Fe(k)  (0,T)  "  Ae(k)  (4) 

From  (1),  (2),  (3)  and  (4)  we  obtain 

G  =  F  (T)  +  F  (T)  +  P.(o,  T,  £, a  )  +  P9(a,  T,£  a) 

°i  °2  1  * 

~Apl  -  AP2  (5) 

The  entropy  production  rate  due  to  propagation  of  the 
CCZ,  as  outlined  above,  is  simply  expressed  as 


T  Si  = 


.  £  + 


Here,  £  and  a  are  the  thermodynamic  fluxes  conjugated  with  the 

“  3  G  ^  3  G 

thermodynamic  forces:  X.  =  \ g  and  X  =  —  .  These  forces 

£  3£  a  9a 

may  be  represented  as  follows 

X£  =  J  +  °£  -  Yeff  *  a  (7) 

and  X  =  M  +  D  -  y  *,#  •  £ 
a  a  'eff 


with  J  =  - 


being  the  well  known  energy  release  rate, 


M  =  - 


is  a  similar  parameter  which  also  has  path-integral 


3a 


representation  for  linear  media.  D„  and  D  are  the  total 

la 

dissipation  energy  rates  associated  with  l  and  a,  respectively. 

An  effective  surface  energy  y  ^  consists  of  the  sum  of  the 
difference  between  the  free  energy  density  of  crazed  and 
uncrazed  materials  (fog  -  fQ^),  crac^  surface  energy 

density  yc  and  the  density  of  the  potential  energy  of  crazed 
material  p^,  i.e., 

Yeff  =  fo2  ~  f o 1  +  Yc  +  P1  (8) 

Generally  y£  is  a  negligible  quantity  in  comparison  to  the  free 
energy  difference  (f  g  Due  to  the  smaYl  stress  in  the 

vicinity  of  the  free  edges  of  a  crack,  p^  is  also  a  small 
quant ity . 

Conventionally,  the  law  of  crack  propagation  is 
established  by  defining  the  relationship  of  the  thermodynamic 
force  ,  represented  by  stress  intensity  factor  K,  or  the  energy 
release  rate  J,  and  critical  surface  energy  density  y,  to  the 
conjugated  flux  described  as  the  crack  growth  rate  l.  Because 
only  one  kinematic  parameter  is  usually  considered,  i.e.,  £, 
the  constitutive  law  has  been  always  approached  by  invoking 
the  first  law  of  thermodynamics  (by  satisfying  the  energy 
balance).  When  two  or  more  kinematic  parameters  are  considered, 
the  energy  balance  alone  would  not  suffice  for  the  deduction  of 
a  constitutive  law.  In  our  case,  where  l  and  a  are  considered, 
an  additional  principle  is  necessary.  Therefore,  a  specific 
version  of  the  second  law  of  thermodynamics,  that  is,  the 
principle  of  minimum  thermodynamic  forces  [7],  is  considered. 


In  view  of  the  second  law  of  thermodynamics,  the  entropy 
production  rate  has  to  be  non-negative.  In  the  simple  case  of 
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a  =  0,  we  can  obtain  from  (6)  and  (7): 

*  • 

T  Si  =  (J  +  -  Yeff  ,ot^  11  -  °  (9) 

Consequently,  the  thermodynamic  forces  associated  with  crack 
* 

growth  rate  (Si  _>  0)  must  also  be  non-negative,  i.e., 

J  +  D£  -  y  .  a  _>  0  (10) 

Experimental  observations  of  fracture  surfaces  [1] 
clearly  show  that  crack  propagates  through  the  root  craze  by 
distinct  jumps  in  a  discontinuous  fashion  (discrete  crack 
advance).  Between  jumps,  damage  within  the  craze  fibrils 
accumulates  until  the  system  reaches  a  state  of  instability 
at  which  another  jump  occurs.  Examination  of  the  parallel 
discontinuous  growth  bands  of  many  polymers  indicates  that  a 
single  jump  approximates  the  length  of  a  root  craze.  In  this 
connection,  it  is  believed  that  a  new  root  craze  evolves 
instantaneously  past  each  jump. 

Since  both  the  energy  release  rate,  J,  and  the  effective 
surface  energy  Yeff  depend  only  on  CCZ  configuration,  therefore, 
they  will  not  change  between  consecutive  crack  jumps  where  the 
system  is  considered  stable.  The  dissipative  energy,  on  the 
other  hand,  accumulates  during  load  excursions.  We  may  there¬ 
fore,  write: 

D£(AN)  =  ¥(£, a, a,  T)  •  AN  (11) 

where  1  is  the  energy  dissipated  per  cycle  and  AN  is  the  number 
of  cycles  per  crack  jump. 

According  to  the  second  law  of  thermodynamics,  crack 
may  grow  only  when  the  condition  (10)  is  satisfied,  i.e., 

J  +  <F  .  AN  -  Yeff  .  a  >  o  (12) 
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This  expression  describes  necessary,  but  insufficient 
conditions  of  instability  for  the  general  case.  Nevertheless, 
under  the  particular  conditions  considered,  where  a  uniform 
stress  is  applied  on  the  boundary  of  the  system,  (12)  is 
necessary  and  sufficient.  When  instability  is  reached,  i.e., 
equality  of  (12)  is  achieved,  CCZ  jumps  into  a  new  configuration, 
at  which  case 


AN 


'eff 


a  -  J 


(13) 


As  we  alluded  to  earlier,  the  length  of  a  crack  jump 
(Ail)  is  considered  equal  to  the  length  of  the  root  craze.  The 
average  speed  of  CCZ  propagation  in  the  direction  of  the  root 
crack  is  therefore  given  by: 


A£  T  •  A£ 


AN 


Yeff  '  a  -  J 


(14) 


In  the  following  section,  we  proceed  to  evaluate  the  independent 
variables  of  this  equation. 


EQUATION  OF  CCZ  PROPAGATION  IN  FATIGUE 
A .  Energy  Dissipation  in  the  Root  Craze 

In  order  to  calculate  the  energy  dissipation  per  cycle 
in  the  root  craze  4*,  we  introduce  few  modifications  to  the 
Dugdale  model  [8]  which  is  illustrated  in  Fig.  2.  In  accordance 
with  experimental  observations  indicating  discontinuous  crack 
growth  [6],  we  suggest  that  the  root  craze  does  not  undergo  any 
significant  propagation  in  the  period  between  two  crack  jumps. 
During  this  period  of  no  growth, the  fibrils  are  assumed  to  be 
under  a  finite  stress,  distributed  uniformly  and  defined  as 
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o*  (t)  #  0.  The  stress  distribution  assumed  (Fig.  2)  provides 
the  following  relationship  between  the  crack  length  (£)  and 
root  craze  (2b)  in  terms  of  external  stress  o^t)  and  o+(T): 

2b  ir'a»(t) 

-  =  2  sin  ( - )  (16) 

l  4  ( t ) 

It  should  be  noted  that  the  formalism  of  the  relationship  (16) 
is  identical  to  that  derived  by  Dugdale  [8].  However,  the 
physical  meaning  of  t)  is  different  from  that  associated 

with  0^  in  Dugdale' s  model  as  explained  later  in  the  discussion. 

The  microscopic  evidence  obtained  from  fatigue  crack 
growth  in  many  polymers  justifies  the  assumption  of  a  constant 
2b/ i  ratio,  within  the  period  between  two  consecutive  jumps. 
Accordingly,  and  from  (16): 

o  (t) 

oo 

-  =  constant  (17) 

For  sinusoidal  loading,  the  time  dependent  stress  may  be 
expressed  as  the  following  function: 


a  (t)  = 


max 


[(1  +  R)  +  (1  -  R)  sin  wt ] 


(18) 


o  . 

where  R  =  - —  is  the  load  ratio  and  w  is  the  load  frequency. 

a 

max 


From  (17)  to  (18),  we  may  write 


a^(t)  = 


*max 


t(l  +  R)  +  (1  -  R)  sin  wt] 


(19) 


Here,  c*max  relates  to  the  critical  craze  initiation  stress 
which  formally  corresponds  to  Dugdale' s  oy.  This  stress  is 
considered  as  the  stress  at  which  a  "new"  root  craze  initiates 


ahead  of  the  crack  at  the  moment  of  crack  jump. 
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During  a  load  cycle,  the  craze  is  believed  to  maintain 
the  same  geometry.  Hence,  the  craze  opening  displacement  is 
described  by  the  following  function: 

<5  ( x ,  t )  =  6*(t)  •  <p(x)  (20) 

Here,  <5*  is  the  opening  displacement  at  the  crack/root  craze 
boundary  (Fig.  2).  The  function  <J>(x)  expresses  the  normalized 
distribution  of  the  craze  opening  displacement  along  its  axis 
and  is  equal  to  unity  at  the  crack/root  craze  boundary.  To 
account  for  the  dissipation  processes  occurring  within  the  craze 
material,  it  is  reasonable  to  introduce  a  phase  lag  X  of  the 
displacement  <S  +  (t)  with  respect  to  o^t).  In  accordance  with 
equation  (19),  we  may  write: 

6*(t)  =  5*max2  [d  +  R)  +  (1  -  R)  sin  (u)t  -  X]  (21) 

From  which  the  rate  of  crack/craze  opening  displacement  is 

1-R 

6.(t)  =  6*  •  -  •  u)  cos  (cot  -  X)  (22) 

*  *max  2 

The  dissipation  energy  per  cycle  *1  of  equation  (11)  can 
be  given  as: 


2tt/o)  2b 

V  =  /  /  o*(t)  ’  S 

o  o 

which  upon  integration  yields: 

6 


*(x,t)  dx  •  dt 


V  61  G*max 


♦max 


(23) 


(24) 


The  term  Bj  of  the  above  expression  is  given  by: 

.2  2b 


(1-R) 


01  = 


sin  X  •  /  4>(x)dx 

o 


(25) 


In  view  of  recent  developments  [9]  in  which  the  Dugdale 
crack  opening  displacement  was  derived  in  terms  of  the  stress 


a 


intensity  factor,  we  may  similarly  write: 


°*max 


^  *max 


max 


E 


(26) 


Finally,  the  energy  dissipation  associated  with  the  root  craze 
per  cycle  can  be  given  by 


K2 

max 

E 


(27) 


B.  Dependency  of  Crack  Jump  on  the  Stress  Intensity  Factor 

We  may  recall  that  the  crack  propagates  through  the  root 

craze  by  distinct  jumps,  each  such  jump  is  =  2b  (the  root 

craze  length).  From  (16),  and  (17),  it  is  obvious  that  tsl  is 

proportional  to  the  crack  length  l  (Fig.  2).  In  fracture 

mechanics  terms,  the  stress  intensity  factor  is  related  to  the 

2  2 

crack  length  through  the  relationship:  K  -  a  tt fc.  V(£/B), 
where  the  function  V  reflects  the  dependence  of  the  stress 
intensity  factor  on  the  geometry  of  the  specimen.  B  is  a 
characteristic  scale  and  is  usually  taken  as  the  width  of  the 
specimen.  Consequently,  the  crack  length  is  given  by 


11  =  B2  K  max 


(28) 


where  1  =  a2  tt  •  V(^/B).  As  £/B<<l,  the  function  V(K/B^  1 

dZt 

2 

and  the  crack  length  is  proportional  to  K  . 

Using  equations  (14),  (25)  and  (26)  and  substituting 


K 


J  = 


max 


,  the  equation  of  CCZ  propagation  may  ultimately 


take  the  form 


E. 


K 


max 


Y 


ef  f  • 


a  -  K 


2 

max 


(29) 


where  3 


St  S0  • 


DISCUSSION 


The  term  E*Y*eff  .a  of  equation  (29)  can  be  considered 

equivalent  to  the  conventional  critical  stress  intensity  factor 
2 

Kc>  Under  conditions  where  a  may  vary  during  crack  growth 
2 

history,  Kc  may  become  a  history-dependent  parameter. 

Three  distinct  crack  growth  regimes  can  be  recognized 
according  to  equation  (29).  The  first  regime  is  defined  by 


K  <<  1  at  which  case  equation  (29)  can  be  rewritten  as: 

K2 

c  4 

d£  $  ^  max  K  2  $  . 

-  =  . .  [l  +  ( — —  —  )  +  ...l'v  2  K 

«  K2  Kc  K 

c  c 


(30) 


Evidently,  this  regime  is  well  approximated  by  Paris  equation 


[11].  The  condition  K  ^  1  defines  the  third  regime  in  which 

K2 


d£  2 

^  vs.  K  exhibits  an  asymptotic  behavior  (Fig.  3).  Regime  II 
can  be  quite  complex  as  shown  by  the  dotted  lines  in  the  figure. 
The  complexity  of  this  regime  is  strongly  dependent  on  changes 
of  a  which  may  be  envisioned  negligible  in  regimes  I  and  III  as 
will  be  explained  below. 

Solution  of  equation  (29)  for  two  different  values  of 
a(a-^  <  otg)  is  represented  in  Fig.  4,  where  the  strong  dependence 
of  the  crack  growth  behavior  on  the  magnitude  of  a  is  shown. 

It  is  therefore  argued  that  the  critical  stress  intenstiy  factor, 
conventionally  thought  to  be  a  material  property,  could  be 
strongly  dependent  on  the  CCZ  width  and  related  changes.  Such 
dependence  has  been  established  theoretically  and  experimentally 
for  other  materials  [10].  In  that  report,  Kc  was  found  to  be 
strongly  dependent  on  the  history  of  loading.  It  should  be 
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emphasized  here  that  the  growth  of  a  represents  only  one  source 
for  changing  Kc>  Th  evolution  of  another  microcrazing  zone  far 
ahead  from  the  CCZ  constitutes  another  influential  source  to 
affect  the  magnitude  of  Kc«  Naturally,  the  model  proposed  here 
can  account  for  such  diverse  phenomena.  During  the  crack 
propagation  history,  a  changes  into  a  higher  value  thereby  causing 
a  reduction  in  the  crack  propagation  rate  or  even  crack 
deceleration.  This  case  is  represented  by  the  dotted  line 
bridging  the  two  constant  a  curves  (Fig.  4).  Such  a  plateau  has 
been  repeatedly  observed  in  several  materials  including  polymers 
under  fatigue  loading  conditions.  An  example  is  shown  in  Fig.  5. 
The  observed  deceleration  in  crack  growth  strongly  indicates  the 
existence  of  a  mechanism  of  energy  absorption  different  from  the 
main  crack  growth.  Elucidation  of  such  mechanism(s)  is  important 
for  understanding  factors  governing  material  toughness. 


APPLICATION  TO  DATA 

As  an  illustrative  test  for  the  model,  some  fatigue  crack 
growth  data  [12]  are  shown  in  Fig.  5  in  the  following  form  of 
equation  (29): 


(31) 


2 

Since  the  width  of  CCZ  was  not  reported,  the  term  K£  has  been 

2 

approximated  as  an  average  of  3.0(MPa)  .m  from  the  recently 


accumulated  data  of  polystyrene  [13].  The  dependency  of  8  on  R 

was  estimated  from  the  experimental  data  of  figure  5.  These 

values  are  tabulated  below: 

R  0  0.2  0.3  0.4  0.6 


6  12  7.5  4.0 


1.56  0.46 


Inspite  of  the  crude  estimation  made  to  evaluate 
equation  (29),  due  to  the  absence  of  precise  measurements 
of  the  related  parameters,  the  model  provides  a  reasonable 
description  of  the  fatigue  crack  propagation.  In  a  future 
publication,  we  report  on  more  accurate  measurements  of  the 
CCZ  zone  in  various  polymers  together  with  further  refinement 
of  the  model. 
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FIG.  1  ( 

FIG.  2 

FIG.  3 

FIG.  4 


FIGURE  CAPTIONS 
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Composite  micrograph  showing  microcrazing 
associated  with  crack  growth  in  polystyrene 
reproduced  from  ref.  [3],  and  (b)  Model  of 
the  crack-crazing  zone  (see  text). 


Crack-root  craze  configuration  according  to 
Dugdale  model.  Note  that  a, 

Dugdale. 


replaces  a  of 

y 


An  illustration  of  three-regime  crack  propagation 
history  according  to  equation  (29). 

Crack  propagation  behavior  at  two  different  values 
of  a.  The  dashes  illustrates  the  transitional 
behavior  when  a  changes  into  a  larger 
configuration. 


FIG.  5 


Fatigue  crack  propagation  data  for  polystyrene 
taken  from  reference  [12]  and  plotted  in  terms 
of  the  proposed  model. 
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